Orientational and induced contributions to the depolarized Rayleigh spectra 
of liquid and supercooled ortho-terphenyl 



O 

o 

(N 



o 



c 

o 
o 



(N 
O 
(N 
O 



c 

o 
o 



13 



S. Mossa 1,2 , G. Ruocco 2 , and M. Sampoli 3 
Center for Polymer studies and Department of Physics, 
Boston University, Boston, Massachusetts 02215 
2 Dipartimento di Fisica, INFM UdR and INFM Center for Statistical Mechanics and Complexity, 
Universita di Roma "La Sapienza" , Piazzale Aldo Moro 2, 1-00185, Roma, Italy 
3 Dipartimento di Energetica and INFM, Universita di Firenze, Via Santa Marta 3 , Firenze, 1-50139, Italy 

(Dated: February 1, 2008) 

The depolarized light scattering spectra of the glass forming liquid ortho-terphenyl have been 
calculated in the low frequency region using molecular dynamics simulation. Realistic system's 
configurations are produced by using a recent flexible molecular model and combined with two 
limiting polarizability schemes, both of them using the dipole-induced-dipole contributions at first 
and second order. The calculated Raman spectral shape are in good agreement with the experimental 
results in a large temperature range. The analysis of the different contributions to the Raman 
spectra emphasizes that the orientational and the collision-induced (translational) terms lie on the 
same time-scale and are of comparable intensity. Moreover, the cross terms are always found to be 
an important contribution to the scattering intensity. 

PACS numbers: PACS numbers: 64.70.Pf, 71.15.Pd, 61. 25. Em, 61.20.p 



INTRODUCTION 

Depolarized light scattering (DLS) spectroscopy has 
been proved to be a valuable tool to explore the dy- 
namical properties of molecular fluids In the case 
of supercooled and glass forming liquids, DLS has been 
used by several researchers with the aim to elucidate the 
dynamical mechanisms underlying the liquid-glass tran- 
sition ||, ||, H ||, 0, |j. However, the uncertainties 
in the scattering mechanisms giving rise to the Raman 
spectra, and the consequent hypotheses that have been 
introduced to interpret the experiments, leave open the 
problem of a reliable analysis. 

In some glass forming liquids, DLS spectra have been 
interpreted as arising primarily from interaction-induced 
mechanisms — namely dipole-induced-dipole (DID) - 
and so they have been related directly to the dynamics 
of the density fluctuations and compared with the out- 
comings of mode-coupling theories (MCT) . Following 
this line, Cummins et al. || have interpreted the spectra 
of salol as arising entirely from the DID mechanism; in a 
more recent work however, the authors assume that 
orientational contributions dominate the DLS spectra of 
salol at least up to 4000 GHz, and the connection with 
the density fluctuations dynamics has to rely on a strong 
coupling between rotational and translational degrees of 
freedom. The situation is similar also for other systems. 
In the case of an other prototype of fragile glass-forming 
liquid, the ortho-terphenyl (OTP), Patkoswski et al. |ll[] 
state that the low frequency part of the spectrum, up to 
about 10 GHz, is due to collective reorientation, while at 
higher frequency the DID contribution is dominating. 

We think that a careful investigation of the scattering 
mechanisms which contribute to DLS spectra is needed 



in order to clarify what type of information is possible 
to extract from the experimental data, and to ascertain 
the connections with the density fluctuation correlators 
treated by MCT. 

In the present work we investigate the low frequency 
DLS spectra of liquid OTP for temperatures above and 
close to the MCT critical temperature T c . We use the 
new realistic intra- and inter-molecular potential model 
set up by some of us very recently 113, na, 14, nsl . Vari- 



ous single particle correlation functions have already been 
studied through classical molecular dynamics (MD) sim- 
ulations, yielding noticeable results. It is now possible to 
combine realistic dynamical configurations with plausible 
polarizability models, and try to discriminate among the 
different hypotheses on the origin of the different contri- 
butions to DLS experimental spectra. 

The paper is organized as follows: Sec. is a brief re- 
call of our potential model and of the methods used to 
produce the dynamical configurations. In Sec. we recall 
the main polarizability models that have been used so far 
to calculate the DLS spectra of simple molecular liquids; 
moreover, we show how they can be used to estimate the 
DLS spectra for our realistic flexible model. The compu- 
tational details are treated in Sec. , while the results and 
a discussion are presented in Sec. . In Sec. we address 
some conclusions. 



THE MOLECULAR MODEL 

In our flexible model |L2| the OTP molecule is con- 
stituted by three rigid hexagons (phenyl rings) of side 
L a = 0.139 nm. Two adjacent vertices of the parent 
(central) ring are bonded to one vertex of the two side 
rings by bonds whose length at equilibrium is Lb = 0.15 
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nm. Each vertex of the hexagons is thought to be oc- 
cupied by a hctious atom of mass Mch = 13 a.m.u. 
representing a carbon-hydrogen (C-H) pair. In the iso- 
lated molecule equilibrium configuration, the two lateral 
rings lie in planes that form an angle of about 54° with 
respect to the central ring's plane. The three rings of 
a given molecule interact among themselves by an in- 
tramolecular potential Vi n t ra ■ This potential is chosen in 
such a way to preserve the molecule from dissociation, 
to give the correct relative equilibrium positions for the 
three phenyl rings, and to represent at best the isolated 
molecule vibrational spectrum. In particular, it is writ- 
ten in the form Vintra = Sfc c fc^fc( r )' wnere each term 
Vfe(r) controls a particular degree of freedom; the actual 
values of the coupling constants Ck have been determined 
in order to have a realistic isolated molecule vibrational 
spectrum. Several internal motions have been taken into 
account, like the stretching along the ring- ring bonds and 
between the side rings, and the rotation of the lateral 
rings along the ring-ring bonds. 

The intermolecular part of the potential, concerning 
the interactions among the rings pertaining to different 
molecules, is described by a site-site pairwise additive 
Lennard Jones (LJ) 6—12 potential, each site correspond- 
ing to one of the six hexagons vertices. The details of the 
intramolecular and intermolecular interaction potentials, 
together with the values of the potential parameters, are 
reported in Ref. |l2| . 

It is worth noting that previous studies of the temper- 
ature dependence of the self diffusion constant and 
of the structural (a) relaxation times jl2| indicate 
that this model is capable to quantitatively reproduce 
the dynamical behavior of the real system, but the ac- 
tual simulated thermodynamical point has to be shifted 
by w 20 K upward. In the following, when we compare 
the simulation results with the experiments, the reported 
MD temperatures are rescaled by such an amount. 

POLARIZ ABILITY MODELS 

In a liquid composed of N optical anisotropic units (a 
unit can be the entire molecule, a group of atoms in the 
molecule, or a single atom inside the molecule), charac- 
terized by a permanent polarizability tensor a, the clas- 
sical low frequency DLS spectrum in the dipole approxi- 
mation is proportional — apart from a trivial frequency 
factor — to the Fourier transform of the time correla- 
tion function of the traceless tensorial part II2 of the 
collective polarizability II. II can be approximated by 
the sum of two terms: the first, II M , is the sum of all 
the permanent polarizabilities dependent on the orienta- 
tional variables Qi of the units; the second, U 1 , is due 
to all the interaction induced contributions, and its lead- 
ing part can be written in term of the dipole propagator 
tensor T^ 2 '(i,j) — V 2 i?~ 1 (where Rij is the distance 



between the i-th and j-th units) : 

n = n M + n J (i) 
u M = J2 WQO} = E( al + pQi) ( 2 ) 

i i 

n 1 = Y^Y^atiTWfrjfa (3) 

i 3 

i j k 
+ ••• 

Here II M is written for units with symmetric top sym- 
metry, a and 7 are the isotropic and anisotropic polariz- 
abilities, and Q a p = (3/2)u Q u ( 3 — (l/2)6 a p, where u is 
a unit vector along the symmetry axis. 

Usually only the first order DID is considered in 
Eq. (||), the second order DID being sometimes com- 
pletely disregarded. The relative importance of the 
higher order DID contributions with respect to the first 
one is dependent on the strength of the permanent polar- 
izability tensor, and on the minimum approach distances 
of the units However, it is well known from the 

study of noble gas systems, that at short distances, when 
high order contributions become important, also other 
multi-pole and quantum electron correlation effects have 
to be considered. These terms are usually found to be 
of opposite sign respect to the DID contribution [|l7], [l8| . 
Empirical diatomic induced polarizabilities have negative 
corrections at short distances for all noble gases [|l9| , and 
only for gaseous mercury a small positive correction to 
the second order DID has been used to fit the experi- 
mental spectra [^0|. From these considerations, in the 
present work we discuss only the contributions from the 
first and second order DID. 

We want to underline that, in the case of OTP, the 
choice of the unit is not so obvious as for monatomic 
fluids or simple molecular fluids such as N%. Thus, the 
subsequent choices are based on the following consider- 
ations. Even for simple molecular systems, like CO2 or 
CS2 fluids, we can consider as unit the entire molecule, 
each atom, or the peripheral (most polarizable) atoms. In 
the case more units are pertaining to the same molecule, 
these units are at short distance each other, so their in- 
duced contributions are tentatively treated in a more re- 
fined scheme. The major attempt in this way has been 
made by Applequist [Mj. His formalism, introduced to 
reproduce the molecular polarizabilities with transfer- 
able atomic polarizabilities, is equivalent to take into ac- 
count the DID at all orders. Indeed the distances are so 
short that the major contributions come from higher or- 
der DID interactions. Irrespective to the relative success 
of the method, various authors |22|, |24| have modi- 
fied the scheme adding mono-pole polarizability, and/or 
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using non isotropic atomic polarizability, and/or sepa- 
rating atoms of the same species on the base of chemical 
bonds. Only a few attempts have been made so far to 
use the Applequist formalism to calculate the experimen- 
tal spectra |l(|, and with doubtful improvements. Thus, 
at present we are not confident in using classical dipolc 
interactions at short distances. 

So far we have not tried to use a complex polarizability 
scheme for our molecular model, and we have chosen to 
adopt two models that can give a rough estimate of a 
lower and upper limit of the induced interactions. In 
our potential model, three six sites rigid rings form the 
molecule and this suggests to consider as units the rings 
themselves (R scheme), or the six sites of each ring (S 
scheme). We have always neglected the effects on the 
polarizability due to the change of one (or two) H-C with 
a C-C bond. In the R scheme, the polarizability tensor is 
taken equal to that of benzene (CqHq). In the S scheme, 
we have simply divided by six the polarizability tensor of 
benzene, and we have neglected the DID contributions 
coming from atoms pertaining to the same ring. The 
induced contributions grow up at least with the sixth 
power of the inverse of the distance, so we expect the R 
scheme to be a lower limit to the induced contributions, 
while the S scheme to be an upper one. Indeed, in the S 
scheme the mean distances are shorter than in R scheme. 
We want to underline that, if in the S scheme the DID 
interactions inside the ring were taken into account, this 
would have reduced the polarizabilities of the six sites. 

Using a notation similar to that of Stassen et al. [ p5| , 
we can work out the collective polarizability correlation 
function for the depolarized scattering: 



C{t) 



1 



(n a (t):n a (o)>, 



(4) 



where V is the volume occupied by the N units. C(t) 
can be divided in three terms, namely the orientational 
term C or (t), the interaction- induced or DID term Ca(t), 
and the cross contribution term C cr (t): 



C(t) = C or (t) + C u (t) + C cr (t), 



(5) 



with 



C or (t) = _ 7 2(Q(0):Q(;)} 

C u (t) = ^(HftOjrllftO)) (6) 

C cr (t) = ^ 7 [<Q(0):n*(*)) + <n£(0):Q(t)>]. 

Here Q = J2i Qi an d au the sums are extended over the 
N units. 

Although the importance of the three terms in Eq. (|^) 
has been found similar in various systems the cross 
term C cr (t) has been neglected, without sound argu- 
ments, in various works on glass- forming liquids. The 



integrated intensities of the orientational contribution, 
i.e., the value at t = Q of the correspondent correlation 
function, can be written as: 



C or (t = 0) = ^ 2 (l + f 2 ), 



(T) 



where gi = 1 + fi is the static angular pair correlation 
factor between our units, defined as lEql 



(8) 



i j^i 



It is well known that, if the units are spherically sym- 
metric (7 = 0), the orientational and cross contributions 
do vanish, as in the case of DLS from monatomic fluids. 
In both R and S schemes, we expect a value of 52 rather 
different from 1. Indeed, also in the R scheme we have 
a strong orientational correlation between different units 
(rings pertaining to the same molecule). This is quite 
different from the case of simple molecular fluids, such as 
N%, where only very small orientational correlations ex- 
ist between different molecules, even in the liquid phase 
near the melting point. If all the rings of OTP molecules 
were in the equilibrium positions, while there was no cor- 
relation between the orientation of different molecules, 
we can calculate the 52 factor for both R and S models. 
Using the equilibrium values of the cosines between the 
normals to the ring planes (ni • ri2 = ni • 113 = 0.588, and 
112 • 113 = 0.673, with the index 1 standing for the central 
ring Q), we have 



92 



1.14, 



g 2 = 6 x g 2 



6.84. 



(9) 



Obviously, the obtained gi factors make the orientational 
contribution integrated intensity — in this idealized con- 
dition — independent of the adopted scheme. Since the 
relation g 2 ' eq = g 2 ' eq '/6 holds always, in the following 
we refer simply to (72 in the R scheme. Further, we can 
evaluate the average value of the analogous of gi for the 
rings pertaining to the same molecules, i.e. 



K? 



3 

2N 



EE< cos2 

i jjii 



'<:) 



) - 1/3), 



(10) 



where the sum is extended over the rings of a molecule, 
9ij is the angle between the normals of the i-th and the 
j-th ring, and () means a configuration or time average. 
The value of K 2 for the isolated molecule at equilibrium 
is always K 2 q ~ 1.14. 



COMPUTATIONAL DETAILS 

We study a system composed by 108 OTP molecules 
(324 rings, 1944 LJ interaction sites) enclosed in a cubic 
box with periodic boundary condition. The MD runs are 
performed in the microcanonical (NVE) ensemble. To 
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FIG. 1: (a) First order DID spectrum calculated considering 
the carbon atoms (model S) at T = 420 K together with the 
resolution curve (see text for details). The insets show the 
data, renormalized to the peak intensity, in semi-logarithmic 
(left panel) and logarithmic scale (right panel), in order to 
stress the power law behavior of the tail, (b) DID spectrum 
calculated considering the phenyl rings centers of mass (model 
R). 

integrate the equations of motion we treats each ring as 
a separate rigid body, identified by the position of its 
center of mass Ri and by its orientation, expressed in 
terms of quaternions qi (^tJ. The standard Verlet leap- 
frog algorithm |27| has been implemented to integrate 
the translational motion while, for the most difficult ori- 
entational part, a refined algorithm has been used |p8| . 
The integration time-step is Ai=2 fs, which gives rise to 
an overall energy conservation better than 0.01 % of the 
kinetic energy. 

We consider two series of runs. For temperature T = 
420 K, to be compared with the experimental results, 
we perform one run 24 ns long with a saving time of 1 
ps; using windows 800 ps wide we achieve a resolution of 
0.04 cm" 1 . For the temperatures T = 280, 290, 300, 310, 
320, 330, 350, 370, 390, 410 and 430 K we run for 640 ps 
with a saving time of 0.04 ps; using windows of 320 ps 
we reach a resolution of 0.1 cm . 

In the R scheme, the components of the polarizability 
tensor (in the molecular fixed frame) are an — 12.31 A 3 , 
a± = 6.35 A 3 , a = l/3(a± + 2a { y) = 10.32 A 3 , 7 = 
a± — ay = —5.96 A 3 , where _L stands for orthogonal to 
the phenyl ring's plane, i.e. parallel to the symmetry axis. 
These values are taken from the polarizability tensor of 
the benzene p9fl. For the S scheme, to each site have 




wave number V [ cm 1 ] 
FIG. 2: (a) As in Fig. 1(a). (b) Second order DID calculated 
for model S. 



been attributed the previous values divided by 6 and the 
same orientation. 
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wave number V [ cm 1 ] 
FIG. 3: (a) Sum of first and second order DID terms for the 
S model, (b) Orientational contribution to the spectrum. 
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FIG. 4: Contributions to the total spectrum for the S model: 
first order DID (circles), orientational contribution (trian- 
gles), and cross-correlation (squares) calculated at T = 420 K; 
the inset shows the same data, renormalized to the respective 
peak intensities, in a double-logarithmic scale. 




wave number V [ cm 1 ] 

FIG. 5: Contributions to the total spectrum for the R model: 
first order DID (circles) , orientational contribution (triangles) 
and cross-correlation (squares) calculated at T — 420 K; the 
inset shows the same data, renormalized to the respective 
peak intensities, in a double-logarithmic scale. 

RESULTS AND DISCUSSION 

In Fig. 1 the DID spectrum at first order is reported 
in the liquid phase at 420 K for both S (Fig. 1(a)) and 
R (Fig. 1(b)) schemes. The dashed lines refer to the 
obtained resolution, which depends on the time length of 
the MD run. In the insets the same data, rescaled to the 
peak intensity, are shown in semi-logarithmic (left) and 
logarithmic (right) scale in order to emphasize the power 
law behavior of the tails. 

As expected, the overall intensity in the S scheme is 
much higher than the one in the R scheme. On the con- 
trary, the shape is very similar, although the relative in- 
tensity of the spectral wings (at 10 cm -1 ) with respect 
to the peak intensity is about three times higher in the R 
scheme. As a consequence the absolute wing intensities 
are comparable. On the overall, the shape is similar to 
the induced contribution from heavy noble gases, i.e., an 
exponential decay of the form exp(— lu/u> ) with oj of the 
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FIG. 6: (a) Total spectra for the S (circles) and R (triangles) 
models, (b) Total spectra for the S (circles) and R (triangles) 
models to geth er with the experimental results (Vh spectra 
from Ref. Q) at T = 444 K (diamonds) in a double loga- 
rithmic scale. A power law of exponent A = —1.5 (dashed 
line) is also plotted as a guide to the eye. 

order of 2 cm -1 . 

Only for the S scheme the second order DID is of 
some importance, so in Fig. 2 the corresponding spec- 
trum (Fig 2(b)) is compared with the first order DID 
(Fig 2(a)). At all frequencies the second order contribu- 
tion is a relatively small fraction of the first one; above 
0.25 cm -1 it is less than about 10%, so the DID's at 
higher order can be safely neglected. 

The sum of the first and second order DID, and the 
orientational spectra in the S scheme are compared in 
Fig. 3(a) and (b) respectively. As we can see, the relative 
shapes and intensities in this case are quite similar (apart 
from the shape of the peak at frequency comparable to 
the spectral resolution, where the simulation statistical 
errors are large and prevent any reliable comparison). 
We want to stress that no time scale separation is de- 
tectable between translational and orientational dynam- 
ics. This is reflected also in the shape and high intensity 
of the cross contribution; at this temperature, indeed, 
the spectra of orientational, DID and cross contributions 
practically superimpose each other, as shown in Fig. 4. 

In the case of the R scheme, the situation is quite differ- 
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FIG. 7: (a) Ratios between the first order DID contribution 
to the S model and the other terms as a function of tem- 
perature: Orientational (circles), DID calculated for model 
R (full squares) and second order DID calculated for model 
S (triangles). The arrows indicate the mean values that are 
respectively 1.22, 0.25, and 0.09. (b) Relative intensity of 
the different contributions with respect to the total intensity 
calculated with both models. The arrows indicate the mean 
values of respectively 0.66, 0.32 0.27, and 0.14. 



cnt. Beyond the smaller relative importance of the DID 
and cross terms, at high frequency (above 5 cm -1 ) the 
spectrum is dominated by the DID contributions, even 
if at 10 cm -1 the ratio of DID, cross and orientational 
contributions is not so high, being 4:2:1 (see Fig. 5). 

Again, even in R scheme our simulations are in poor 
agreement with the assumptions made by Patkowski et 
al. ||l"l"f about the spectral separation between orienta- 
tional and DID contributions. Indeed, it is evident from 
Fig. 5 that no time scale separation can be effective and 
the cross term cannot be neglected. 

The resulting — total — simulated DLS spectra for 
the two schemes are reported in Fig. 6(a). In Fig. 6(b) 
we show, in double logarithmic scale, a comparison of the 
two line shapes with the experimental spectrum at 444 
K | p0| . The agreement among the sets of data is quite 
good, considering the approximations we have made in 
the dynamical and polarizability model. This shape com- 
parison gives some preference to the S scheme, but it is 
hard to really discriminate between these two limiting 
polarizability models on this basis. R and S DLS spectra 
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FIG. 8: Static angular pair correlation factor (72 (T) (Eq. (8)) 
as a function of temperature. In the inset, K2IT) (Eq. (10)) 
is shown as a function of temperature. 



differ essentially in intensity rather then in shape, and 
some help would come from an absolute intensity com- 
parison; unfortunately, the OTP experimental absolute 
intensity has not been calculated so far. Anyway, the 
large uncertainties usually present in both experimental 
and simulated determinations of DLS absolute spectral 
intensities, could not allow this kind of comparison. 

In Fig. 7 we plot the temperature dependence of the 
relative integrated intensities of the various spectral con- 
tributions. We see that the relative intensities do not 
show any significant trend with temperature, apart from 
a tendency of the orientational part to increase at the 
lowest investigated temperatures. The increase has to be 
attributed to a change in the relative orientation of the 
rings, i.e., to an increasing of 172- This is clearly seen in 
Fig. 8 where gi (main panel) and Ki Inset) are plotted 
as a function of temperature. It is interesting to note 
that also K2 increases on decreasing temperature, being 
less than the reference value K^ q — 1.14 at high tem- 
peratures and more than that value at low temperatures. 
On the basis of the value of the angles between the sym- 
metric top axis of the rings, K2 increases fast if there is a 
spread of ni -ri2 and ni 113, but decreases with the spread 
of ri2 • 113. The distribution of angle cosines is reported 
in Fig. 6 of Ref. ]l2] ]. At high temperatures, the large 
spread of 112 • 113 is, anyway, able to decrease K2 under 
, but at low temperature the situation is reversed by 
the spread of ni • n 2 and ni • n 3 . 

We expect a similar situation to hold for the DLS of 
a large class of molecular glass forming liquids: the DLS 
would be a mixture of orientational, DID and cross con- 
tributions in the entire low frequency range, with no sig- 
nificant time scale separation. Further, the better agree- 
ment with the experimental results of the S scheme in the 
case of OTP, underlines the importance to take into ac- 
count the internal degrees of freedom to obtain a realistic 
description of the DLS of glass forming liquids consisting 
of large flexible molecules. 
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CONCLUSIONS 

In this paper we have studied, by means of molecu- 
lar dynamics simulations, the orientational and induced 
contributions to the low frequency depolarized Rayleigh 
spectra of supercooled ortho-terphenyl. We have used 
a realistic flexible intramolecular model recently intro- 
duced by some of us, taking into account the most im- 
portant internal degrees of freedom of the OTP molecule. 
Two polarizability models have been introduced, each of 
them considering a different scattering unit. In the ring's 
model R, we have assigned to each phenyl ring the polar- 
izability of the benzene. In the site's model S, we have 
assigned to the sites of each ring the same polarizability 
divided by six. This procedure has allowed us to study 
the single contributions to the DLS spectra in both cases. 
Our main findings are the following: i) Although in the 
two schemes the intensities of the DID contributions are 
very different, their overall shape is very similar; ii) The 
second order DID is at all frequencies already a relatively 
small fraction of the first order contributions, so that all 
the ligher induced terms can be safely neglected: in) For 



[4 
[5 
[6 

ir 

[8 
[9 
[10 

[11 

[12 
[13 
[14 
_[15_ 



the 5 model first order DID and orientational contribu- 



tions are very similar in both shape and intensity; iv) In 
both models the cross correlation between induced and 
orientational terms cannot be neglected. This fact is in 
striking contrast with the analysis of the experimental 
data of Ref. [[fl| and support the conclusion that: v) No 
time scale separation is present between DID and orienta- 
tional contribution. In other words we expect that, sim- 
ilarly to the present case, for a broad class of molecular 
glass forming liquids the DLS spectra would be a super- 
position of DID, orientational and mixed contributions, 
and none of them should be disregarded. Finally: vi) 
the better agreement of the S scheme calculated spectra 
with the experimental results underlines the importance 
of taking into account the molecular internal degrees of 
freedom in order to obtain a realistic description of com- 
plex molecular liquids. 
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